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ON  THE  STATISTICAL  THEORY  OF  TURBULENCE* 
By  W.  Heisenberg 


The  interpretation  of  turbulence  presented  in  the  preceding  paper 
by  V.  Weizsaecker  is  treated  mathematically  with  the  aid  of  the  custom- 
ary method  of  Fourier  analysis.   The  spectrum  of  the  turbulent  motion 
is  derived  to  the  smallest  wave  lengths,  that  is,  into  the  laminar 
region;  the  mean  pressure  fluctuations  and  the  correlation  functions  are 
calculated.   Finally,  an  attempt  is  made  to  derive  the  constant  which  is 
chairacteristic  for  the  energy  dissipation  in  the  statistical  turbulent 
motion  from  the  hydrodynamic  equations. 

In  the  statistical  theory  of  turbulence  developed  by  G.  I.  Taylorl 
and  V.  KarmdJi,2  the  irregular  turbulent  motion  of  a  fluid  is  described 
by  several  characteristic  functions  between  which  simple  mathematical 
relations  exist:   the  "spectral"  distribution  of  energy  to  waves  of  dif- 
ferent wave  length,  the  correlations  between  the  velocities  at  points 
along  a  prescribed  displacement  in  space  or  time,  axid  the  like.   The 
reports  of  G.  I.  Taylor  contain  detailed  empirical  and  theoretical  data 
on  these  functions. 

In  his  preceding  report,  v.  Weizsaecker?  derived  the  most  important 
of  these  functions,  the  spectral  distribution  of  energy  for  the  limiting 

*"Zur  statistischen  Theorie  der  Turbulenz."   Zeitschrift  ftir  Physik, 
vol.  124,  19ij-8,  pp.  628-657. 

iTaylor,  G.  I.:   Proc.  Roy.  Soc.  A  I5I,  ^21  (1935);  156,  50?  (1956); 
l64,  15  (1958);  i6if,  476  (1958). 

2Karman,  Th.  v.:   Journ.  Aero.  Sci.  h,    151  (1937). 

3c.  F.  V.  Weizsaecker:   ZS  f.  Phys.,  being  published.   V.  Weizsaecker 's 
paper  and  the  present  treatise  have  been  written  in  close  collaboration 
during  the  time  of  our  stay  in  England  in  1945.   Only  after  the  articles 
had  been  finished,  Mr.  G.  I.  Taylor  kindly  told  us  (spring  19^6)  that 
essential  ideas  in  these  articles  had  been  found  and  published  already 
by  Koljnogorof f :   Compt.  Rend.  Acad.  Sc.  USSR  30,  301  (19^1);  52,  I6  (19^-1); 
and  Onsager  (Phys.  Revue  68,  286,  19^5).   Compare  a  report  by  G.  K. 
Batchelor  at  the  VI.  Internat.  Kongress  F.  angew.  Mechanik,  Paris  19^6. 
Approximately  at  the  same  time  we  learned,  furthermore,  about  a  paper 
by  Prandtl  and  Wieghardt  which  contains  similar  concepts  and  has  mean- 
while been  published  in  the  GSttingen  Academy  reports  (Nachrichten  der 
Akademie  der  Wissenschaften  in  G8ttingen,  Math. -Phys ikal  KLasse  of  the 
year  19^5,  p.  6).   The  present  paper  may  therefore  be  regarded  only  as 
a  supplement  and  completion  of  these  earlier  investigations. 
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case  of  large  Reynolds  numbers  on  the  basis  of  similitude  considerations. 
The  following  sections  of  this  paper  first  will  translate  v.  Weizsaecker ' s 
considerations  into  the  accustomed  terms  of  Fourier  analysis  and,  with 
the  aid  of  this  translation,  study  the  discontinuity  of  the  spectrum  for 
high  frequencies  due  to  molecular  viscosity.   Then  conclusions  will  be 
drawn  for  the  correlation  functions  and  the  pressure  fluctuations,  and 
finally,  a  derivation  of  the  fundamental  constant  of  the  energy  dissipa- 
tion will  be  attempted. 


1.  Representation  of  v.  Weizsaecker ' s  Considerations 

in  the  Terms  of  Fourier  Components 

At  sufficiently  high  Reynolds  numbers  the  energy  dissipation  for 
the  turbulent  motion  takes  place  in  such  a  manner  that  the  large  turbu- 
lence elements  lose  energy  due  to  the  fact  that,  for  them,  the  energy 
and  momentum  transfer  by  small  turbulence  elements  has  the  effect  of  an 
additional  viscosity  (cf.  for  instance  Prandtl^). 

Under  steady-state  conditions,  energy  thus  is  continuously  trans- 
ferred from  larger  to  smaller  tiirbulence  elements,  with  the  spectral 
region  of  a  certain  wave  length  always  receiving  from  larger  waves  as 
much  energy  as  it  gives  off  to  smaller  waves.   For  maintenance  of  this 
equilibrixim  a  certain  energy  distribution  is  necessary  which,  when  the 
molecular  friction  is  neglected,  is  represented  according  to  v.  Weizsaecker 
by  the  law 


p^  =  p^  =  p  r  F(k)dk  (1) 

F(k)  ~  k-5/5 

=  — -     signifies  the  wave  mamber,     Vp,  =    1/ v       or  measure  of 

^    wave  length  '-'    V 

the  mean  velocity. 5 j  This  spectnim  F(k)  ~  k-5/3  is  bounded  on  two  sides. 
For  small  wave  nurabers,  that  is,  large  wave  lengths,  for  one  reason  or 
another,  it  will  no  longer  be  possible  to  regard  the  flow  as  isotropically 
turbulent,  for  the  largest  t\irbulence  elements  are  governed  by  the  geom- 
etry of  the  devices  which  generate  the  turbulence.   This  end  of  the  spec- 
trum for  small  k  therefore  cannot  in  any  case  be  the  subject  of  a  purely 
statistical  theory.   For  large  k,  in  contrast,  the  spectrum  is  bounded 

^L.  Prandtl,  StrCmungslehre  (Flow  theory).   Vieweg,  Braunschweig, 
3rd  edition  19^2,  p.  105  ff. 

50ur  Vq  differs  from  the  Vq  in  v.  Weizsaecker ' s  report  by  a 

numerical  factor  of  the  order  of  magnitude  of  1. 
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by  the  molecular  viscosity.   For  large  k,  finally,  the  molecular  vis- 
cosity will  become  larger  than  the  apparent  turbulent  viscosity,  and 
the  spectrum  then  will  drop  off  very  rapidly. 

For  the  calculations,  we  shall  use  the  following  notation:   the 
velocity  v  in  a  normalized  volume  V  is  to  be  expanded  into  a  Fourier 
series 


V  -  N  I^^'^L   =  g  n,, 


Ij  "=='  • 


}    nx,   ny,   nz  integersj     (2) 


Therein  Vjj.  =  v^-k  ^^^  ^^^   number  of  the  "natural  vibrations"  between 
k  and  k  +  Ak  is  given"  by  ^^ — ^^,   Then  one  obtains 

(at)5 


o        o   ^ —  — k 


=  i  /lf:rk2dk— L_ 


^ 


2  - 


F(k)dk 


(3) 


thus 


F(k)  =  (2rt)"^^ 


^k 


(^) 


From  div  v  =  0  there  follows 


(%^ 


(5) 


Let  us  call  the  coefficient  of  the  moleculai-  viscosity,   ^l.   The 
mean  energy  loss  then  is,  because  of  div  v  =  0  with  the  assumption 
that  the  bounding  surfaces  are  at  rest: 


°This  method  is  somewhat  less  apparent,  but  mathematically  more 
convenient  than  the  customary  expansion  with  respect  to  sin  an.d  cos. 
It  amounts  formally  to  the  limiting  condition  that  y, Sv/Bx,  ...  at  a 
bounding  surface  of  the  volume  V  are  to  have  the  same  values  as  at  the 
opposite  bounding  surface.  ' 


S  =  ti(rotv)^ 


=  [J^ 


z 


M 


=   n  /  F(k)2k2(ik: 
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(6) 
(7) 


If  the  spectrum  obeys  in  a  large  region  the  law  F(k)  ~  k~5/5,  the 
total  energy  is  determined  by  the  largest  tiirbulence  elements.  We  may 

as  Slime  for  instance  that  the  law  k~5/3  is  valid  down  to  a  minimum  wave 
number  kQ;  for  smaller  k  we  shall  assume  F/,  \  =  0.   Then 


Vp,2  =  2  /   F(k)dk  =  2 


Thus 


and 


ko  k 


c 
573 


dk  =  5CkQ' 


■2/5 


(8) 


F(k)  = 


c=  :^  ^=/5 


vo^lEo!^ 


,     ,   f  or  k  >  kp, 
5   v5/5       -  0 


(9) 


vo2ko2/3(2n)3 
6nVk^V3 


(10) 


V.  Weizsaecker  considers  the  energy  loss  Sjj.  of  that  portion  of 
the  total  spectrum  the  wave  numbers  of  which  lie  below  k.   For  these 
turbulence  elements,  the  turbulence  elements  of  smaller  wave  length 
(<23t/k)  have  the  effect  of  an  additional  viscosity.   One  may  therefore 
generally  write 


Sk  =  (^^  +  ^k)  /   F(k')  2k'2dk 


(11) 


where 


T\,       is  to  designate  the  additional  tiorbulent  viscosity;  it  i^ 

produced  by  the  ciimulative  action  of  all  turbulence  elements  with  wave 
lengths   <2jt/k.   With  respect  to  dimensions,  t]t^     is  according  to 

Prandtl  the  product  of  density,  mixing  length,  and  velocity  where  the 
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mixing  length  will  be  coinpcirable  to  the  diameter  of  the  turbulence  ele- 
ments in  question  whereas  the  velocity  of  the  turbulence  elements  is 

given  for  instance  by  VQ(kQ/k)-^/5.  With  reference  to  v.  Weizsaecker 's 

report  one  will  therefore  put 

Tlj,  =  .p  /   dk'  /^i^  (12) 


(  k;   is  a  numerical  factor) 

The  expression  below  the  integration  sign  is  essentially  determined  by 
dimensional  considerations;  but  one  could  of  course  imagine  that,  for 
instance,  the  waves  k'   in  the  proximity  of  k  enter  into  the  integral 
with  somewhat  different  weights  than  the  waves  with  large  k  values; 
that  is,  the  integrand  could  depend,  in  addition,  on  the  dimensionless 
number  k'/k.   However,  because  of  the  homogeneous  form  of  the  spectrum 

F(k)  ~  k~5/5^  one  may  include  all  these  uncertainties  in  the  numerical 
factor  K   and  give  to  the  integral  arbitrarily  the  exact  form  (12). 

This  method  is  unobjectionable  in  the  region  of  thfe  k"5/^_law  but 
becomes  inaccurate  at  the  ends  of  the  region  where  the  geometry  or  the 
molecular  friction  modifies  the  spectrum.   But  even  at  the  latter  limit 
(12)  will  still  be  a  good  approximation  which  at  least  qualitatively 
correctly  represents  the  effect  of  friction. 

The  constant  k   in  eq^uation  (l2)  must  be  exactly  determined  by 
the  hydrodynamic  equations;  it  has  the  same  numerical  value  in  all  cases 
where  one  may  speak  of  statistical  isotropic  turbulence,  and  does  not 
depend  in  any  way  on  the  geometry  of  the  flow.   The  theoretical  deter- 
mination of  this  important  number  will  be  attempted  in  section  5-   There, 
it  will  also  be  shown  that  the  turbulent  energy  dissipation  actually  may 
be  written  as  a  double  integral  of  the  type 

/^k  O  CO 

dk'.../   dk''...,  with  the  integrand  signifying  the  energy 
0  ^k 

transferred  from  k'   into  k' '   per  unit  time  (equation  (89)).   This 
integral  is  more  complicated  than  the  simplifying  expression  (ij)  which 
results  from  (ll)  and  (12);  however,  for  the  following  considerations  (ll) 
and  (12)  may  be  regarded  as  sufficient  approximations.   For  S-^     one  thus 
obtains 


(13) 


The  decisive  step  of  v.  Weizsaecker 's  consideration  is  the  state- 
ment that  this  expression  for  k  »  kQ  must  be  independent  of  k: 
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Sjj^  =  S  =  const  (for  k  »  Rq)  {ih) 

because  the  total  energy  lies,  almost  entirely,  in  the  long-wave  region 
of  the  spectrum,  and  the  energy  "transport"  thus  must  become  independent 
of  k. 

Equation  (ij)  may  therefore  be  interpreted  as  the  determining  equa- 
tion of  the  turbulent  spectrum  F(k)  which  must  yield  for  the  region  of 

large  Reynolds  niombers  the  k-5/5-law,  and  for  still  larger  k  values 
the  fading  of  the  spectrum  due  to  the  molecular  viscosity. 


2.  Shape  of  the  Spectrum  in  the  Region  of  the 

Smallest  Turbulence  Elements 

We  put  first  n/p  =  v  and  differentiate  (l3)  with  respect  to  k. 
Then  there  results 


^+  dk'>/£^JF(k)  k2  =  l/FCkl  /   F(k')k'2dk'       (15) 


Then  we  define  new  variables  x  and  w  by  the  equations 

X  =  Ig  A,  F(k)  =  F(ko)e-^;  w  =  w(x)  (l6) 


Therewith  (l5)  is  transformed  into 


(17) 


The  constant  ^  1^^  =   — :r^7Z~     ^'^^'    equation  (9))  is  essentially  the 


'fo    ^vq 

reciprocal  Reynolds  number  of  the  total  flow  and  therefore  always  very 
small;  if  the  Reynolds  number  itself  were  small,  the  flow  could  not  be 
turbulent  at  all.   By  repeated  differentiation  there  originates  from  (l?) 
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In  this  equation^  one  may  approximately  evaluate  the  integral 
poo    _H±2i 
/   dxe      by  expanding  w(x)   in  the  proximity  of  x:  ^(^1 )  = 

X 

w(x)  +  fX2_-x\   -^  +  .  .  .   and  breaking  off  after  the  second  term.   Since 

the  exponential  function  rapidly  decays,  one  thus  obtains  a  good  approx- 
imation 

w+x 
000    _  w+x      — — 

/   dxe "  2  ~  2  -S ( 19 ) 

^x  1  +  |2 

dx 

By  substitution  into  (18)  there  finally  results 

7  _  ^V— 1-^  e"^  +    1   U  2  (20) 

dx/ 

One  recognizes  from  (20)  immediately  the  variation  of  the  spectrum.  For 
not  too  large  x  and  w  the  first  term  in  the  sum  may  be  neglected  and 
one  obtains 

(^  -  i)  =  <^  -  s>      '-  i = f'  ^(^) = 'o{i^f'     (21) 

as  must  be  the  case  according  to  v.  Weizsaecker ' s  theory.   For  large  x 
and  w,  in  contrast,  the  first  term  predominates;  therefore  one  must 
then  have 

^  =  7   F  ~  const  k"'''  (22) 

dx 

In  the  region  of  the  smallest  turbulence  elements  the  spectrum  therefore 
decays  very  rapidly,  namely  with  the  seventh  power  of  the  wave  number. 

Only  in  the  transitional  region  from  (2l)  to  (22)  are  numerical 
calculations  necessary  in  order  to  determine  the  solution  of  (20).   Since 
one  may  put  for  smaller  x,  that  is,  for  the  region 

1  «  X  «J  Ig  2L\ 

4    V  yko 

(23) 

5 

w  =  -f  X 

5 
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(therewith  not  only  (l8),  but  also  (17)  then  is  satisfied  with  sufficient 
accuracy),  one  may  calculate,  progressing  from  point  to  point,   dw/dx 
according  to  (20)  from  w,  and  therewith  derive  w  for  higher  x.   It 
is  siifficient  to  perform  the  numerical  calculation  for  a  particular 


large  value  of  the  constant,  for  instajice  —  ,  —^  =  a. 


For  another 
value  b  one  may  then  obtain  w  by  a  simple  similitude  transformation: 


w^(x)  =w^fx  +|ls|)  -ilg| 


{2k) 


as  one  recognizes  by  substitution  into  (20)  and  (23). 

Figure  1  shows  the  result  of  the  numerical  calculation  for 

2k  /Fn 

—  /— ^  =  about  1000.   The  numerical  integration  shows  that  for  this 

value  in  the  region 

X  >  5 

w(x)  ~  7x  -  21.85 

More  generally,  one  obtains  therefore  in  this  region  of  the  k~ '  -  law: 


(25) 


that  is. 


w(x)=7x^5.0-4lg|j|^ 


3  k  /     . 
F(k)  =  o.oksG  ^(^)  fe 


(26) 


(27) 


A  serviceable  interpolation  formula  which  is  correct  in  the  two 
limiting  cases  and  does  not  result  in  any  large  errors  in  the  transitional 
region,  either,  reads: 


F(k)  =  ?o\^j 


U5/3 


^M& 


.\«/' 


J 


(28) 


If  one  defines  I^  =  n/kg  as  the  "diameter  of  the  largest  turbulence 
elements"  and  introduces  as  the  Reynolds  number  of  the  total  flow 


Rr 


PVqIo 


(29) 
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one  obtains  according  to  (9)^  (2?)^  and  (29) 


,6ko(RoK)5/^  (50) 


kg  =  o.ie 

One  may  denote  L3  =  it /kg  as  the  "diameter  of  the  smallest  turbulence 
elements"  and  obtains 

Ls  =  6.25Io(Ro^)"^/^  (51) 

By  (9)^  (28),  and  (30)  the  form  of  the  spectrum  is  determined  for 
the  entire  k-region.   For  the  actual  flows,  of  course,  the  shape  of 
the  spectrum  will  be  different  for  small  k-values   (k~kQ)   since  there 
the  geometry  of  the  tests  plays  a  role,  for  instance,  the  shape  of  the 
grids  by  means  of  which  the  turbulence  is  produced.   In  order  to  be  able 
to  sensibly  carry  out  the  comparison  with  experiments,  one  will  then 
introduce  a  quantity  kQ  in  such  a  manner  that,  for  instance  in  the 

domain  of  the  k"^/^-law  (thus  for  kg  »  k  »  kg),  the  formila 

becomes  correct.   The  quantity  kg  thus  defined  then  does  not  give  any 

direct  statement  regarding  the  variation  of  the  spectrum  for  the  smallest 

k-values.   Generally,  though,  the  spectrum  will  greatly  deviate  from  the 

-5/3 
k    -law  in  the  region  k  ~kQ. 

For  k  »  kg,  too,  the  spectrum  will  not  unlimitedly  retain  the  form 

k~ '  .   The  well-known  investigations  of  Burgers '  mak:e  it  very  probable  thai 
at  sufficiently  small  Reynolds  numbers  finally  no  turbulent  motions  what- 

_7 
ever  exist.   On  the  other  hand,  the  k  '-law  shows  such  rapid  decay  that 

the  region  k  »  kg   is  practically  insignificant.   A  somewhat  larger 

error  will  arise,  particularly  in  the  transitional  region,  due  to  the 

inaccuracy  of  equation  (ij)  itself j  but  it  is  probably  not  worth  while 

to  apply  aljready  at  this  point  the  much  more  complicated  equation  of 

section  5  "to  the  problem  stated  here.   The  correct  equations  would  at 

any  rate  lead  to  somewhat  different  numerical  factors  in  (27),  (30),  and 

(31). 

'('Burgers,  J.  M.  :   Verb.  d.  Kgl.  Nied.  Akad.  d  Wiss.  I7,  Nr.  2,  1 
(1939);  18,  Nr.  1,  1  (19^0). 
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For  the  comparison  with  experiment  one  needs  the  energy  distribution 
with  wave  rnmber  for  a  certain  direction,  for  instance  with  k^,  since 

the  spectra  have  been  measiored  experimentally  by  Simmonso  and  Dryden9 
by  means  of  the  fluctuations  of  the  velocity  with  time  in  an  airstream 
which  is  guided  past  the  measuring  point  at  a  constant  velocity  U  which 
is  large  relative  to  Vq.     Also,  this  spectrum  has  a  different  form 

according  to  whether  the  Fourier  expansion  of  v^  or  vy  is  concerned. 
Experimentally,  first  the  spectrum  for  v^  is  reqiiired;  however,  we 
shall  also  derive  the  spectrum  for  v   since  it  will  be  necessary 
later  on,  in  the  calculation  of  the  correlation  functions.   Since 
according  to  equation  (5)   (Xkis)  -  ^^    '-"^®  obtains 

The  spectrum  of  v^  -^"^  •'Sc  ^^^^^^  ^^  shall  designate  by  Fx(k^) 
becomes  therefore 


^x(kx) 


:r^if-§->' 


1  /    dk  (vii  _  V  ^lF(k)  (5^) 


'4,  ?^'  -"'' 


In  a  similar  manner,  the  spectrum  for  Vy  becomes 


F 

y 


wv:r^if-^K> 


From  (3^)  and  (35)  there  follows  for 
(a)  kQ  «  k  «  kg: 

(36) 

8  ' — ' — 

Simmons,  U.  Salter:   Proc .  Roy.  Soc .  A  I65,  73,  (1938). 

-^Dryden,  Schubauer,  Mock  u.  Skramstad:   Nation.  Adv.  Comm.  Aero., 
Nr.  581,  1938;  Dryden,  H.  L. :   Proc.  V.  Intern.  Congr.  f.  Applied.  Mech. 
Cambridge,  Mass.,  p.  362,  I938. 
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("b)  k  »  kc 


F(k)  =  F, 


0 


1^5/3^316/5 


kT 


,{^)   -  gg  ko5/\,W5,^-T 


•y 


M' 


2F(- 
63" 


1^ 


5/3v  16/3,-7 


(37) 


As  a  serviceable  interpolation  formula  (which,  however,  is  somewhat 
less  accurate  than  (28)  in  the  transitional  region),  one  may  again  put 


with 


■\^}  -  note 


\^xj 


1  + 


k.  \8/3 


.k.^y 


-2 


(38) 


XX 


=  0.6k^  kc 


(39) 


and 


•w  -  S^fef 


55     Ux/ 


1  + 


.^s^y 


8/5 


(i+O) 


with 


k/""  =  0.793  k. 


(^1) 


Before  the  comparison  with  experiment  is  carried  out  in  detail,  we  want 
to  raise  the  question  at  what  critical  Reynolds  numbers  the  transition 

from  the  ^'^'-^-lav   to  the  k~'^-law  takes  place,  that  is,  -  if  one  wants 
to  express  it  in  this  manner  -  the  transition  from  the  really  turbulent 
motion  proper  to  the  laminar  motion.   One  may  regard  as  the  critical 
Reynolds  nimber  for  this,  perhaps,  the  expression 


Rs  = 


wherein,  according  to  v.  Weizsaecker,   Vg 

and  (31)  one  then  obtains 

R= 


10.2 


(^2) 


/k  \i/3 
vo(^)    .   From  (29),  (30), 


{^3) 
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The  numerical  value  of  k.  will  be  discussed  later  on.   At  any  rate,  the 

traxLsition  therefore  takes  place  at  a  certain  numerical  value  of  the 

Reynolds  number,  as  was  to  be  expected  from  general  similitude  considera- 
tions. 

In  figure  2,  the  measiorements  of  the  spectrum  F^(k^)  by  Simmons^Q 

are  compeared  with  the  theory.   The  measurements  in  question  ajre  intensity 
measiirements  on  an  airstream  which  flows  past  the  measuring  point  at  the 
velocities  U  =  h^6   cm/sec  (O),  6o8  cm/sec  (X)^  IO60  cm/sec  (^]) ,    and 
which  has  been  made  tiirbulent  by  a  grid  of  7-6  cm  mesh  width;  the  meas- 
urement was  made  2.1  m  behind  the  grid.   The  measured  points  of  Simmons 
are  plotted  individually  only  in  the  right  part  of  the  figure,  in  the 
left  half  the  approximate  scatter  of  the  measured  points  is  indicated 

by  a  vertical  line.   The  abscissa  is  k  in  cm"-'-  (in  logarithmlcal 
scale),  the  ordinate  F^fk^^),  likewise  logarithmically,  In  arbitrary 

units.   If  one  assumes  U/vq  to  have  the  same  value  in  all  three  meas- 
uring series,  which  is  confirmed  by  other  measurements  by  Taylor,  one 
obtains,  in  the  case  of  a  suitable  selection  of  this  ratio,  the  three 
curves  plotted  in  the  figure.   Qualitatively,  the  Simmons  data  are  well 
represented  by  the  curves,  particularly  also  the  divergence  of  the  three 
test  series  in  the  short-wave  part  of  the  spectrum.   In  details,  however, 
there  exist  considerable  discrepancies;  one  recognizes  from  the  figure 

that  the  range  of  validity  of  the  k'-^'-^-law  is  here  so  small  that  a 
reliable  check  is  not  possible.   The  reason  for  this  is  the  smallness 

of  the  Reynolds  number  Rq.   For  k  =  1  cm"   the  diameter  of  the  turbu- 
lence elements  is  3  cm,  thus  about  half  as  large  as  the  mesh  width  of 
the  grid)  in  this  region,  the  turbulence  is  not  yet  fully  Isotropic, 

therefore  the  k"-'/^  -law  cannot  yet  be  valid.   However,  already  for 

k  =  4  cm"-'-  the  influence  of  the  molecular  viscosity  becomes  noticeable, 
and  the  intensity  drops  off  markedly.   The  related  measurements  by  Dryden 
quoted  before  which  extend  over  a  large  spectral  region  have  been  made 

at  Reynolds  numbers  so  small  that  the  validity  of  the  k~5/5_]_aw  can 
hardly  be  checked.   Therefore  it  would  be  desirable  that  similar  meas- 
urements be  carried  out  at  very  much  larger  Reynolds  numbers .   For  the 
ratio  U/vq,  one  obtains  from  the  adjustment  of  the  theoretical  curves 

to  the  measiirlng  points  of  U/vq  =  53i^   if  one  identifies  I^  with  the 

mesh  width  of  the  grid.   This  value  agrees  well  with  measurements!!  of 
this  ratio  in  similar  tests  if  one  assumes  k      to  be  about  0.5' 

l^Simmons,  U.  Salter:   Proc.  Roy.  Soc.  A  165,  73  (1938). 
llCf.  G.  I.  Taylor,  Proc.  Royal  Soc.  A,  l6k,    k-Q6,    (1938). 
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Another  and  probably  more  accurate  determination  of  k      is  obtained 
from  the  damping  with  time  of  the  turbulence,  already  theoretically  treated 
completely  by  Taylor  (cited  before).   For  the  total  energy  loss  per  cm5 
and  second:   S,  there  results  from  (9)^  (l3)^  and  (l^) 

S  =  p.  J^  Vq^o  (i+^) 

For  the  damping  with  time  of  Vq  there  must  thus  apply 


^^4)...M,^^  (,5) 


with  the  solution-'-^ 


^0 


(0) 


V, 


(t)  = -^ (46) 

1  +  -^K  koVo(0)t 


Taylor  who  essentially  derived  this  eq_uation  reports  on  measurements 
by  Simmons  in  which  U/vQ(t)  was  ascertained  as  a  function  of  t  =  x/u 

(x  equals  distance  of  the  measuring  point  from  the  grid).   From  (46) 

follows 

U       U     \/3  X      U     „  ^Q   X  ,|,„v 

— 7— r  =  — 7— r  +  -^  ^^  —  =  — 7— r  +  0.68k  —         (47) 
Vo(t)   Vo(0)    8     I^   Vq(0)         Lq 

If  one  puts  Vq  =  u' y  3  (u'  =  V  v-j^   according  to  Taylor)  and  identifies 
I^  with  the  mesh  width,  there  follows  from  Taylor's  measurements  k  =  O.85, 
from  the  corresponding  measurements  by  Dryden  a  somewhat  smaller  value. 
However,  because  of  the  uncertainty  regarding  the  value  to  be  inserted 
for  Lq,  this  determination  is  probably  still  uncertain  by  about  50  per- 
cent. 


3.  The  Correlation  Functions 

Taylor  and  von  Karman  ( cited  before )  studied  the  correlations  which 

exist  between  the  velocities  at  two  points  at  a  given  distance.   The  two 
__ 

Footnote  at  the  time  of  proof  correction:   For  this  solution, 

kQ  =  const,   is  presupposed  which  certainly  is  not  the  case  for  larger 

times.   The  problem  of  damping  is  investigated  more  closely  in  a  paper 

of  the  author  about  to  be  published  (Proc.  Roy.  Soc.  A.). 
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correlation  functions  R]_(x)  and  R2(x)  which  therein  play  the  main 
role,  are  defined  as 


Ri(x).M^^4^^)     B,(.) .  M^iVyi'^) 


(M) 


2  2 

with  the  point  Pp  displaced  with  respect  to  the  point  P-]_  by  the 
distance  x  in  the  x-direction. 

These  fimctions  are,  according  to  Taylor,  in  a  simple  relationship 
with  the  spectra: 


Rl(x)  =  ^-JT^ 

/    dkxFxC^x) 


/   dk;xFy(kx)cos  k^x 
R2(x)  =  -^ 

^  0 
With  the  aid  of  equations  (j^J-)  and  (35)^  (^9)  is  transformed  into 

p  00 

5  /    dkF(k)(sin  k:c  -  kx  cos  kx)k~5x"5 
Ri(x)  =  -^ 


/   dkF(k) 

p  /    dkF(k)(k2x2sin  k  x  +  kx  cos  kx  -  sin  kx)k"5x~5 
R2(x)=_^ p (50) 

/   dkP(k) 
^0 

From  these  expressions  one  recognizes  immediately  the  correctness  of 
von  Karman's  relationship 

I^2  =  %  +  |h^  (51) 


NACA  TM  U3I  15 

The  formulas  (50)  may  be  approximately  evaluated  In  the  two  limiting 

cases  X  «  l/kg  and  l/kg  «  x  «  l/kQ.   When  x  «  :r~,    it  is  advisable 

^s 
to  expand  the  integrands  with  respect  to  powers  of  x.   The  first  terms 
of  the  expansion  then  lead  to  the  quantity 

pOO 

/   dkF(k)k2 

k2=10 (52) 

/   dJiF(k) 
-'0 

which  may  easily  be   calculated  from  (15)^    (^^)^    an<3.  (29): 


k2   =  ^KRoko2  (55) 

Thus   one  obtains: 


For     x  «  l/ks  (5^) 

R2(x)  =  1  -  2M  +  .  .  . 
Taylor  defined  a  length  A  by  the  equation 

A^  =  ^  (55) 

and  designated  it  as  a  measure  for  the  magnitude  of  the  smallest  turbu- 
lence elements.   According  to  (53);   ^  becomes 

1 

X  =  2.71  -^  =   0.i^3^Ls(RoO^  (56) 

It  must  be  emphasized  that  A  is  not  identical  with  the  quantity 
Lg   (equation  (51))  which  we  have  denoted  as  "diameter  of  the  smallest 

turbulence  elements"  and  that  A  also  depends  on  I^  and  Vq  in  a 
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manner  different  from  that  of  Lg-   A  comparison  of  (31)  and  (56)  shows 
that  for  sufficiently  large  Reynolds  numbers  Lg  becomes  «A. 

In  the  opposite  limiting  case  l/kg  «  x  «  I/^q  °°-®  obtains 
from  (50) 

Rl(x)  =  1  -  0.6U3(kox)2/5  +  .  .  . 


RgCx)  =  1  -  0.858(kox)2/5  + 


(57) 


Here  these  first  terms  of  the  expansion  do  not  depend  on  the  special  form 
of  the  spectrum  in  the  proximity  of  k^j  only  for  x  ~  l/kQ  does  the  form 

of  the  spectrum  in  the  proximity  of  kg  becomes  important;  there,  however, 

the  problem  may  no  longer  be  treated  with  purely  statistical  methods. 
The  formulas  (5^)  snd   (57)  thus  give  essentially  a  complete  description 
of  the  correlations  in  so  far  as  they  may  be  regarded  as  a  consequence 
of  statistical  isotropic  turbulence. 

The  formulas  (5^)  and.  (57)  also  show  clearly  that  the  correlation 
function  does  not  have  the  same  form  in  all  flows  but  that  rather,  in  the 
case  of  variations  of  the  parameters,  the  inner  and  outer  parts  of  the 
function  undergo  different  similitude  transformations.  This  point  has 
been  stressed  particularly  by  TaylorlJ  in  contrast  to  a  different  con- 
jecture of  von  Karman  (cited  before). 

For  comparison  with  experiment,  the  measurements  of  R2_(x)  and  R2(x) 
made  by  Simmons  have  been  plotted  (circles  and  dots,  respectively)  in 
figure  3;  furthermore,  the  theoretical  curves  calculated  according  to 
the  exact  formula  (50)  are  shown.   Here  again  Lq  has  been  identified 

with  the  mesh  width  of  the  grid  and  A  has  been  calculated  from  the 
spectrum  for  U  =  IO60  cm/sec.   The  experimental  values  agree,  in  fact, 
with  the  theoretical  ones  very  exactly  at  the  smaller  values  of  x, 
actually  more  exactly  than  could  have  been  expected  in  view  of  the  uncer- 
tainty of  I<3.   Beginning  from  xkQ  ~  1  the  deviation  of  the  experimental 
points  from  the  theoretical  curves  becomes  noticeable  which  was  to  be 
expected  from  the  derivation.   The  variation  for  larger  x-values  depends 
on  the  behavior  of  the  spectrum  in  the  proximity  of  kQ  which  cannot, 
in  principle,  be  represented  by  our  formulas.   But  even  for  larger  x-values 
the  deviations  from  the  theoretical  curves  remain  small. 

15g.  I.  Taylor,  Jour.  Aero.  Sci.  k,    No.  8,  51I,  1957. 
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h.    The  Pressure  Fluctuations 

While  studying  the  diffusion  in  a  turbulent  airstream,  Taylor  (cited 
before)  has  derived  a  relationship  between  the  correlation  function  for 
the  diffusion  and  the  root  mean  square  value  of  the  pressure  gradient. 
We  shall  investigate  therefore  also  the  root  mean  sq_uare  values  of  the 
pressure  fluctuations  from  the  viewpoint  of  the  theory  here  described. 


With  reference  to  equation  (2)  one  may  expand  the  pressure  into  a 
Fourier  series 


=   Z-^-e^-T.  .   =  T1  * 


P^e 


^P-k  =  Pk 


(58) 


and  the  fundamental  hydrodynamic  equation 


V  =  -  (vv)v  -  i  VP  +  ^  Av 


is  transformed  into 


^  = 


Z_  i(zk'4ik-k'  -  T  %  "  ^  ^^i 


(59) 


Because  of      Ivj^k]   =  0     it  follows 


ik=  -Z  H 


-         k'  ~ 


k 


-k-^'    "  ^  ^k-k'^^ 


p        -k 


Pk  =   -  ;3Z_(Yk'^)(:^k-k'^k) 


(60) 


For  the  root  mean  square  values  of  pressure  and  pressure  gradient  there 
results 


?  =  I 


k 


Pk 


grad^p 


J_ 


Pk'= 


(61) 


We  are  interested,  first,  in  this  latter  mean  value: 
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grad2  P  =    X5I    Z]  ^  (Yk'^)(Xk-k'^)(i^k"^)(^k-k"t)  (62) 

k    k'      k"   ^  ~ 

The  superscript  bar  indicating  approximation  to  the  mean  signifies 
here  simply  approximation  to  the  mean  in  time.   If  one  wants  to  calculate 
mean  values  of  the  type  (62),  one  must  take  as  a  "basis  some  kind  of 
"ass\imption  of  disorder"  regarding  the  turbulent  motion.   One  may  start 
from  the  fact  that  the  amplitudes  y^     in  the  course  of  time  fluctuate 

by  a  value  given  by  (lO)  or  (28),  respectively,  so  that  the  time  average 

V,  V    is  given  simply  by  the  spectrum  (28).   The  phases  of  the  v,  , 

however,  will  in  the  co\irse  of  time  pass  through  all  possible  values;  all 
values  of  the  phase  will  occur,  on  the  average,  with  the  same  frequency. 
If  one  could  regard  the  phases  pertaining  to  different  wave  numbers  as 
completely  independent  statistically,  there  would,  in  taking  the  mean  of 


such  products  of  four  factors  v,  v,  v,  y,  ,  be  left  only  the  terms  in 
which  every  two  wave  numbers  are  equal  and  opposite;  thus,  terms  of  the 


type  V,  V  ,  y^  v   ,  and  these  mean  values  could  be  replaced  by  the  pro- 
ducts  of  the  mean  values  of  the  squares  of  the  amplitudes: 


Actually,  however,  statistical  correlations  will  exist  between  the  phases 
pertaining  to  different  wave  numbers  since  the  waves  mutually  influence 
one  another.   In  section  5  we  shall  attempt  to  estimate  such  correlations 
in  a  simple  case.   In  spite  of  the  presence  of  the  correlations,  however, 
we  are  inclined  to  assume  that  in  a  sum  of  the  kind  (62)  the  terms  of  the 
type  (65)  make  the  largest  contribution  since  their  mean  value  is  differ- 
ent from  zero  even  in  the  first  approximation,  without  any  assumption 
regarding  the  behavior  of  the  waves  with  respect  to  time  whereas  the 
other  mean  values  attain  a  value  different  from  zero  only  because  of  the 
finer  fluctuations  of  various  waves.   We  believe  therefore  that  one  obtains 
a  serviceable  approximation  if  one  takes  only  the  terms  of  the  type  (63) 
in  (62)  into  consideration.   Then  there  results 


(\'^")(^-K-^")(\"^')(^-k'^') 


Srad^p  =  2p2  /_ l^r IXJl^ 1\^ L}^^ L  (g)^) 

k'k' '  (k'    -  k")2 
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In  taking  the  mean  with  respect  to  the  directions  of  the  v-^,  -  one  assumes 
again  all  directions  perpendicular  to  k  for  Vv  to  be  equally  probable  - 
it  is  expedient  to  use  the  relation 


Yk 


(vka)(Y_kb)  =  — ^ 


ab 


(a^)(:bk) 

k2 


(65) 


Thus  one  obtains 


jrad^)  =  — 


k'k' 


Kk' 


k'^k' '2  -  (k'k' ')2]' 
(k'  -  k")2 


{SG) 


If  one  sets   (k'k' ')  =  k'k' '^,  one  can  perform  the  integration  with  respect 
to  ^   and  finds 


7^  =  fZ^    K'2|  K..2|k'k"J|n-) 

d      3^,j^,,   1-11-    1         \J£   / 


(67) 


where 


.(s)=,i  =-i^ 


1  +  ^  s2  4-  ^  s^  -  s6  +  (^  -  ^^)   Ig 
3       5  2s      *= 


1  +  s 


|1  -  s| 


For  0  <  s  «  1  there  applies  approximately 

^-(3).8£fi-5  32^s^' 
15  \    7     21 


(68) 


(69) 


If  one  transforms  the  sums  into  integrals  and  substitutes  the  spec- 
trum (28)  into  (67),  there  follows  finally 


k 


grad2p=  |p2voV      % 


dk' 


'(^ 


2  r 


5    ^0    ;,,,3 


8" 

2 

8 

2- 

/k'\5 

-1        /k'  ' 

\^ 

1  +     ^^\ 

1  +     - — 

■ 

\_        Vksj  J 

[             Us 

/ 

(70) 
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One  recognizes  from  (TO)  that  the  integrals  converge  at  small  k-values, 
and  that  one  may  therefore  perform  the  integration  from  k  =  0  without  a 

considerable  error.   This  shows  that  grad  p  is  actually  determined  "by  the 
hehavior  of  the  spectrum  at  large  k,  that  is,  by  the  smallest  turbulence 
elements.   Had  we  calculated  p2^  we  would  have  found,  on  the  contrary,  that 
the  integral  diverges  at  small  values  of  k,  th\is  that  its  value  is  deter- 
mined entirely  by  the  largest  turbulence  elements.   Therefore,  the  value  of 
p2  cannot  at  all  be  calculated  according  to  the  method  used  here;  for, 
first,  in  the  case  of  small  k,  the  spectrum  has  a  form  dependent  on  the 
geometry;  and  second,  it  would  surely  by  quite  unjustified  to  consider,  for 
the  largest  turbulence  elements,  only  the  mean  values  of  the  type  (63)  "^ince 
the  geometry  certainly  impresses  definite  phase  relations  upon  the  system  of 
the  largest  eddies . 

Equation  (70)  now  becomes 
gradp-^v^ko   k3    j^  ^^/^J^  ;?75^i,^8/3|2Y,,,8/3j2     ^'^^ 


The  double  integral  at  the  right  was  estimated,  according  to  a  gra- 
phical method,  to  be  O.763;  thus  there  follows  finally  (cf.  (30)) 


grad^p  =  0 .  17p2vQ\QV5kg2/3  (^2) 

=  0.05p2vo\o2y^  (^5) 


Taylor  (cited  before)  had  expressed  the  conjecture  that  grad2p 
should  have  the  same  order  of  magnitude  as   p2vq2(Sv/Sx)2  -thus  the  order 
of  magnitude   p^vg  kQ^/^kg^/^.   One  recognizes  now  from  (72)  that  grad2p 
miist  be  considerably  smaller,  the  more  so,  the  larger  the  ratio  kg/kg. 
The  length  Ar)  defined  by  Taylor: 


1/  -  -^^i 
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that  is^ 

^  k 

grad2p  =  I  p2  ^  (75) 

5    A^2 


must,  therefore,  for  large  Reynolds  numbers,  become  considerably  larger 
than  the  length  >v  of  equation  (55).   From  (75),  (73),  and  (56)  follows 


A^  =  A0.1+2^RoK  (76) 

It  is  true  that  this  result  does  not  agree  with  the  experimental 
findings.   Taylor  indicates,  for  a  test  which  Simmons  had  performed 
adjoining  similar  experiments  by  Schubauer^^,   A^iA  =  0.5  where  one  must 

ass\mie  approximately  yRo'^  ~  3-9  according  to  the  test  conditions. 

Thus  one  must  raise  the  question  whether  the  result  (72)  has  perhaps  been 
falsified  by  the  fact  that  only  the  terms  of  the  type  (63)  were  taken  into 
consideration  in  taking  the  mean.   However,  one  can  easily  see  that  this 
may  perhaps  affect  the  numerical  factor  in  (72)  but  that  the  dependence  of 
kg  and  kg,  that  is,  the  dependence  of  Ati  on  (Rq*^  )   ^^  ^^  ^°   ^^y  con- 
nected with  this  approximation.   For  already  equation  (62)  shows  that  on 
the  right  side  the  normalization  factor  Vq  k^  '    must  appear,  because 

of  equation  (10)  .  After  the  mean  has  been  taken,  this  factor  is  supple- 
mented by  a  factor  of  the  dimension  k2/3  which  obviously  can  be  at  most 
of  the  order  kg2/3  j  it  must  be,  because  the  pertaining  integral  with 
respect  to  k  would  diverge  like  k^/J  if  the  decay  of  the  spectrum 
would  not  set  in  for  k  ~  kg  with  k"7.   There  would  remain  the  possi- 
bility that  only  the  numerical  factor  in  (72)  has  been  estimated  as  too 
low  due  to  the  consideration  solely  of  the  terms  of  the  type  (63).   But 
it  is  hard  to  imagine  that  the  correct  expression  would  increase  by  more 
than  a  tenfold  -  which  would  be  necessary  for  interpretation  of  the 
experiments . 

Perhaps  the  contradiction  may  be  cleared  up  in  the  following  manner: 

The  main  contribution  to  grad^p  stems  from  wave  numbers  of  the  order  kg, 

thus,  from  turbulence  elements  whose  diameter  measures  a  few  millimeters. 
In  Simmons'  test,  the  airstream  is  heated  by  a  hot  wire  of  20  cm  length 
stretched  across  the  wind  tunnel,  and  then  the  distribution  of  the  heated 
air  is  measured  at  a  certain  distance  behind  the  hot  wire .   Precisely  the 
smaller  distances  (5  to  I5  cm)  are  decisive  for  the  determination  of  'K   . 

^^Schubauer:  Rep.  Nat.  Adv.  Coram.  Aero.  Nachr.  Nr.  524,  1955- 
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It  suggests  itself  to  assume  that  the  hot  wire  itself  produces  in 
the  airstream  a  small  vortex  street  and  additional  turbulence,  with  the 
turbulence  elements  probably  having  a  length  of  a  few  millimeters  -  that 
is,  the  wire  increases  precisely  the  intensity  of  the  turbulence  in  the 

spectral  region  which  exerts  the  strongest  influence  on  grad^p.   The 
additional  turbulence  in  the  immediate  proximity  of  the  wire  is  probably 
much  more  intense  than  the  original  turbulence  of  the  same  wave  range . 
However,  this  additional  turbulence  is  rapidly  damped,  of  course,  and  it 
is  siirely  difficult  to  estimate  whether  this  turbulence  alone  can  explain 
the  discrepancy  between  (76)  and  the  empirical  A^n-value. 

5 .  Energy  Dissipation  for  Normal  Isotropic  Turbulence 

The  investigations  of  the  preceding  section  are  already  closely  con- 
nected with  a  basic  problem  of  the  statistical  turbulence  theory:   namely, 
with  the  determination  of  the  energy  dissipation  in  the  case  of  normal 
energy  distribution,  that  is,  the  determination  of  the  constant  k   in 
equation  (12).   In  this  problem  the  molecular  friction  may  be  neglected 
entirely.   The  fundamental  hydrodynamlcal  equations  could  therefore  be 
presupposed  in  the  form 


V  =  -(vV)v  -  i  Vp     Vv 
P 


Furthermore,  a  suitably  singled-out  partial  volume  of  the  fluid  is  to  be 
selected  as  normalization  volume  which  under  certain  circumstances  is 
moved  simultaneously  with  the  fluid,  corresponding,  for  instance,  to  the 
mean  value  of  the  velocity  with  respect  to  the  volume.  We  assume  therefore 
that  the  volume  moves  with  the  velocity  u_.   Then  (60)  is  transformed  into 


V, 


=  -iXI  (Xk'k) 


k  ^ 

-k-k-  ";^^^-k-^) 


i(^)v^       (77) 


For  the  calculation  of  the  energy  dissipation,  one  has  to  ascertain 


how  the  intensity 


v^2   of  a  certain  natural  vibration  (or  perhaps 


better:   the  sum  of  such  squares  of  amplitudes  with  respect  to  a  small 

k+Ak 


spectral  region  Ak:/   Vu   I  varies  in  the  course  of  time.   One 
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recognizes  from  (77)  that  one  needs  for  this  purpose  time  averages  with 
respect  to  products  of  the  type 


%1  ^   %3  (t8) 


wherein  k^z  =  -k-]_  -  kg.   Because  of  the  statistically  uniform  distribution 

of  the  phases,  these  mean  values  would  disappear  if  there  would  not  exist 
statistical  correlations  between  the  phases  pertaining  to  different  k 
which  stem  from  the  mutual  influence  of  the  various  waves,  as  has  been 
explained  already  in  section  k.      In  order  to  ascertain  these  correlations, 
one  must  somehow  express  in  the  equations  the  influence  exerted  upon  a 
wave  with  given  k  by  waves  with  a  different  k;  one  may  do  this  for 
instance  by  representing  one  of  the  three  amplitudes  in  (78)  as  a  time 
integral  over  vj^  and  expressing  v^     in  turn  by  a  sum  over  two  other  vj^ 
according  to  (77)'   Then  one  obtains  products  of  four  amplitudes  v,   each 

of  which,  however,  must  partly  be  taken  at  different  times.  For  such  pro- 
ducts the  considerations  of  section  h   are  valid  according  to  which  one 
obtains  a  first  approximation  by  taking  only  products  of  the  type  (63) 
into  consideration.   Of  course,  one  could  continue  the  procedure  in  prin- 
ciple and  attempt  to  calculate  the  other  mean  values  of  quadruple  products 
by  tracing  them  back  to  six-fold  ones,  etc.   However,  such  calculations 
would  probably  become  much  too  complicated;  the  higher  terms  probably  also 
would  make  a  lesser  contribution,  and  we  shall  thus  be  content  with  the 
first  step. 

For  these  calculations,  one  will  obviously  need  mean  values  of  the 
type 


and  we  define  therefore 


Vk(t)v_j^(t  +  t) 


k+Ak 


Rk(t,T)  =  ^ ^^^     (79) 
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The  simmiation  over  a  small  spectral  region  Z^  has  been  included 
into  the  definition  of  Ri^(t,T)   so  that  the  magnitude  of  the  normaliza- 
tion volume  does  not  directly  enter  into  Rj^(t,T)   and  that  the  mean  is 
taken  equally  over  all  directions  of  k.   Evidently  the  spectral  region  Ak 
must  be  selected  wide  enough  that  many  natural  vibrations  of  the  normaliza- 
tion volume  still  may  be  accommodated  in  it  (that  is,   k^Z^V  »  1),  yet 
very  small  compared  to  k  itself.   These  requirements  are  for  the  turbu- 
lence elements  of  the  order  of  magnitude  V  themselves  no  longer  compatible 
but  for  those,  the  statistical  methods  cannot  be  applied  anyway.   The  whole 
procedure  thus  can  be  carried  out  only  if  it  is  found  that  the  large  turbu- 
lence elements  practically  do  not  any  more  contribute  to  the  mean  values  to 
be  investigated. 

In  order  to  obtain  from  the  hydrodynamical  equations  information 
about  the  quantities  Rj^(t,T),  it  suggests  itself  to  examine  the  following 

expression: 


"^i.H) .  w„.^^^„^-? 


%(*  ^  i)^-kt  -  i)  =  -i ^'v%'^'V  Wl'^-k  j  ^  ^^^^^^  ^-/     ^^°) 


In  this  expression  one  can  replace  v-^t   by  a  time  integral  with  respect 
to  vj^t  : 


t+l-T' 

dT'^,^         +Vi,,(-T)      (81) 


If  T  is  selected  sufficiently  large,  the  correlation  between 
v_j^i(t)   and  Vj^t(T)  will  disappear;  it  is  therefore  expedient  to  per- 
form, after  substitution  into  (80),  the  limiting  process  T — >«=.   If,  in 
addition,  one  takes  the  mean  with  respect  to  the  .lirections  -  thus  elimi- 
nating the  term  with  u  -  one  obtains  from  (80)  and  (8l): 
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(k  k' 


t+^-T' 


-  u^'-k"   ^ 


k" 


Xk-k'^ 


-k" 


(82) 


If  one  furthermore,  as  in  section  h,    takes  into  consideration  the  terms 
of  the  type  (63)  only,  also  replaces  k  -  k'  by  k',  and  integrates  with 
respect  to  the  cosine  of  the  angle  between  k  and  this  vector,  there 
follows : 


k+Ak. 


hi^  -  i)^-k(^ 


k+Z& 
k 


dk'k' 


(2«) 


5  16 


dT' 


vS        vg,    R, 


o'O 


^ 


k(t     -^,T     -     T') 


Rl^.^t  +  ^^,T')k-5k-3(k2    -  k'2)    2kk'(k^  +  k'^    -   |k2k'2)    - 


(k2  +  k'2)(k2    -  k'2)lg  kiJll 


k  -  k' 


(83) 


This  equation  presents  the  possibility  of  expressing  the  differential 

dRv(t,T)       dR^(t,T) 

quotients  with  respect  to  time  — and  — by  the  Rv   them- 

dt  dT 

selves;  when  the  Rj^  are  known,  one  can,  moreover,  calculate  the  energy 

dissipation  from  (83),  putting  t  =  0.   For  this  purpose  we  shall  assume 
that  the  entire  turbulence  phenomenon  is  either  stationary  or  is  damped 

noticeably 


The 


very  slowly  so  that  the  times  during  which  the  intensity  Vj^'^ 

varies  are  very  long  compared  with  the  fluctuation  period  of  v.  . 

notation  Y.]^    \    I'epresents,  therefore,  the  mean  value  over  a  time  which  is 

certainly  much  longer  than  the  fluctuation  period  but  is  very  much  shorter 
than  the  damping  time . 


dR^(t,T) 
The  equation  for  — gives  a  measure  for  the  fluctuations  of 


dt 


the  quantity  Ri^(t,T)   as  a  function  of  t  about  its  time  mean; 


RkC-^)  =  Rk(t.T) 


(84) 
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One  may  surmise  that  these  fluctuations  are  small  in  the  region  of  small  t 
which  is  determined  by  the  small  turbulence  elements,  and  that  they  increase 
with  growing  t;  this  question  will  be  further  investigated  later  on. 

Before  carrying  out  the  further  calculations  one  has  to  determine  how 
the  partial  volume  V  and  its  velocity  u  are  to  be  chosen.   One  could 
first  try  to  put  u  =  0  and  to  identify  V  with  the  total  volume.   How- 
ever, one  would  obtain  an  erroneous  picture  of  the  actual  conditions:   The 

decrease  with  time  of  the  correlation  function  Rj^(t,T)  =  Rj^(t)   as  a 
function  of  t  is  determined  in  this  coordinate  system  by  the  largest 
turbulence  elements,  and  is  therefore  very  rapid.   One  can  show  that  the 
correlation  function  in  this  coordinate  system  is  given  with  sufficient 
approximation  by 


Rk(T)  =1 :=r^ 


kVQT 


^2^^2  pf^^Ol 
n/T2 


dxe-'^ 


2 


The  calculations  which  have  led  to  this  expression  need  not  be  discussed 
in  more  detail  since  the  expression  is  not  further  used  later.   The 
physical  interpretation  of  the  expression  is  given  by  the  following  con- 
sideration:  The  function  Ric(t)   in  it  decreases  after  a  time  of  the 

order  -^^;  that  is  the  time  during  which,  for  instance,  precisely  an 

kVQ 

eddy  of  the  wave  length  — ,  due  to  the  high  velocity  in  the  largest  tur- 

bulence  elements,  passes  by  the  point  of  observation.   The  fact  that  the 
correlation  function  decreases  after  that  time  signifies  therefore  simply 
that  the  velocity  in  the  largest  turbulence  elements  is  of  the  order  of 
magnitude  vq,  but,  statistically,  fluctuates  about  values  of  this  order. 
This  phenomenon  is  not  connected  with  the  disintegration  of  the  eddy  of 

the  wave  length  ^.   If  the  F(k)  ~  k"5/5_iaw  is  valid,  it  is  rather  to 

be  expected,  according  to  the  similitude  considerations  of  v.  Weizsaecker, 

that  the  disintegration  of  the  eddy  takes  place  only  after  time  intervals 

-1  -2/5   -I/5 
of  the  order  Sitvq  k    ^q  '  .   On  the  other  hand,  however,  the  energy 

dissipation  is  connected  with  the  disintegration  of  the  eddies,  not  with 
the  motion  on  a  large  scale.   If  one  wants  to  describe  the  disintegration 
of  the  eddies  in  equations,  one  must  move  the  coordinate  system  at  the 
same  time.   One  then  must  make  the  linear  dimensions  of  the  partial  vol- 
ume V  somewhat,  but  not  very  much  larger  than  ^  and  move  it  simul- 
taneously  in  proportion  to  the  mean  velocity  within  it.  We  shall  assume 
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experimentally  that  one  can  give  for  every  k  a  volume  V  corresponding 

to  it  in  such  a  manner  that  Vk^  becomes  independent  of  k  and  that  for 
the  thus  selected,  simultaneously  moved  volume  elements,  the  correlation 

function  R(t,T)   is  a  universal  function  of  the  variables  VQk  '^kQ' ^t 

and  VQk  /^kQ/^T  as  is  to  be  expected  according  to  v.  Weizsaecker's 

similitude  consideration.   We  shall  show  that  the  relations  for  Rj^(t,T) 

arising  from  (85)  actually  can  be  satisfied  by  this  assumption  if  the  Fj^^ 

are  distributed  according  to  the  k'-^/^-law. 

If  one  puts  Rj^(t,T)  =  g(^,Ti)  wherein 


t  -  IOj,2/3,  1/5.    ^  _  IOv2/3a/5_    k-  _  ^ 


(85) 


there  follows  from  (83) 


dT]  16 


(iyf(y) 


dTi' 


i^ 


IL. 


;({^  +  ^^)y^/5,V: 


■2/3 


g(c  -  ^,T,  .  v)g((^  -  ^)y2/3,Vy^/5J 


d^        "8 


rioo  Poo 

dyf(y)  dTi' 

0  Jo 


g(^  -  ^,^'  -  n 


gf(^  +  ^4-^)y2/5,Vy2/5) 


n+vx  2/3    .2/3 


(5-^,n.v)g(^-^)y^/^Vy 


(86) 


where 


f(y)  = 


1  -y2 


Mh 


2y(l  -  |y2  +  yM 


(1   -  y2)(l   -  y^)lj 


1  +  y 


1  -y 


(87) 
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These  equations  actually  do  no  longer  contain  the  constants  Kq,   Vq. 

The  reason  is  that  the  integral  with  respect  to  y  converges  for  small  as 
well  as  for  large  y;  f (y)  disappears  sufficiently  for  y  =  0  as  well  as 
for  y — >oo.  Therefore  one  may  take  the  integral  over  k',  Instead  of  from 
]fiQ,    simply  from  0,  without  considerable  error;  moreover,  the  convergence 

of  the  integral  for  large  values  of  k  shows  that  the  molecular  friction 
actually  is  of  no  importance  in  this  problem;  the  behavior  of  the  spectrum 
in  the  region  of  the  smallest  turbulence  elements  is  unimportant  for  the 
correlation  functions  R(t,T)   and  the  energy  dissipation  for  medium 
k-values . 

Before  attempting  the  numerical  solution  of  (86),  we  shall  use  the 
equation  (85)  for  calculating  the  energy  dissipation  in  the  approximation 
here  aspired  to.   For  this  purpose  we  put  t  =  0  and  integrate  the  equa- 
tion (83)  with  respect  to  k  between  two  arbitrary  limits  K]_  and  K25 


"IK2 


4 
2 

2    n 

lt(2«)5  Jki 

>                 /Too 

dk            dk' 

^  Jo     ^' 

i 

2 
Xk' 

(k^    -   K.^) 

2kk'fk^  +  k'**  -  |k^l 

[^   dT.R,(t-^,v)R,,(t-^,x.) 


,2 


)-  (k2.k'2)(k2.,.2)\g. 


.2\^.J^  -   k'l 


k  +  k' 


(88) 


The  integrand  on  the  right  side  is  an  antlsymmetrical  function  in  k 
and  k'.   If  one  calculates  the  variation  with  time  of  the  total  energy, 
that  is,  if  one  puts  K]_  =  0,  K2  =  °°,    there  results  therefore  zero,  as 

far  as  the  integral  on  the  right  converges  at  all.   That  is,  the  total 
energy  is  constant  in  time;  this  is  a  necessary  requirement  since  the 
molecular  friction  has  not  been  taken  into  consideration.   However,  if 
one  considers  the  variation  with  time  of  the  energy  which  is  contained  In 
the  part  of  the  spectrum  lying  between  Kj_  and  K2J  the  integral  on  the 

right  side  may  be  transformed  in  the  following  manner  (we  shall  call  the 
antlsymmetrical  integrand  simply  J) : 


i+rtk^dk 


dk 


dk'J 


dk'J 


dk 


dk'(-j) 


(89) 
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In  the  first  of  the  two  integrals  on  the  right  J  is  always  posi- 
tive; in  the  second  (-J)  is  always  positive.   From  this  notation  there 
follows  that  the  first  integral  may  be  interpreted  as  the  energy  which, 
per  unit  time,  flows  from  smaller  wave  numbers  (k'  <  K^)  into  the  region 

between  K.-^     and  K2;  the  second  integral  as  the  energy  which  flows  toward 

larger  wave  numbers  (k'  >  K2).   If  one  puts,  in  particular,  K-,  =  0  and 

Kp  »  k-Q,    the  second  integral  represents  the  entire  energy  dissipation; 

for  the  normal  spectrum  [F(k)  ~  k    J  it  must  prove  to  be  independent  of 
Kg.   Thus  one  obtains  from  (85),  (86),  and  (89)  for  the  energy  dissipa- 
tion the  expression 


Kp 


s  =  ^I'-o 


J^>^;"  a^(.) /;  a,.(,  -!„)«((,  -  |;^/3,,^/3) 


=  |-v5ho  r  dy(-lgy)f(y)  j  "  drig  (^  -  \,^B^    -  f)y^h,^Y^/^  (90) 

This  expression  is  actually  independent  of  Kg  as  it  must  be.   Since 
the  entire  energy  dissipation  according  to  equation  (44)  is  also  given  by 

pK;-^v^k_,  there  follows 
o  u  (J 


cirig(^  -|.n)s((^  -  |)y2/5,T,y2/3^       (91) 


From  this  equation  k;   can  be  calculated  numerically  if  the  fionction 
g(^,Ti)   is  known. 

We  now  turn  to  the  treatment  of  the  equation  system  (86).   This  sys- 
tem represents  a  considerable  simplification  compared  to  the  initial 
equation  (77)  in.  so  far  as  it  does  not  any  more  contain  any  dimensional 

-5/3 
quantities,  and  is  already  derived  from  the  equilibrium  spectrum  k  '  . 

On  the  other  hand,  (86)  also  still  contains  statements  regarding  the 
fluctuations  of  the  g(^,'n)   as  a  fimction  of  the   ^  and  is,  for  that 
reason,  doubtlessly  too  complicated  to  permit  rigorous  solutions.   One 
could  attempt  to  completely  neglect  the  fluctuations  in  a  first  approxi- 
mation, and  to  calculate  with  the  mean  values.   Unfortunately,  however. 
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it  turns  out  that  the  contribution  of  the  fluctuations  in  certain  regions 
is  large.   One  recognizes  this  from  the  second  equation  (86). 

For,  if  one  puts 


g(^;Ti)  =  g(Tl)  +  Ag(^,Ti] 


(92) 


where  g(T\)      signifies  the  mean  value  over   ^  of  g(^,i 


g(-n)  =  g(C.Ti) 


(95) 


there  follows  from  the  time  mean  of  the  second  equation  (86) 


I        ciyf(y)  dV|_g(V    -    Ti)   +   g(V    +   Ti)]g(Ti'y2/3\ 


_      5 


dyf(y)  dT)' 


^K 


2  ' 


V    -   ri]Agf(^   +  IL^jy2/5,Vy2/3V 


Ag(^    -  ^,V    .    ri)Ag[(c    -  l|_Zi:)y2/5,,.y2/5) 


(9M 


This  relationship  shows  that  the  fluctuations  Z^(^,ri)   cannot 
always  be  small.   It  is  true  that  the  left  side  of  (9k)    vanishes  for 
Tj  =  0;  this  follows  from  the  relation 


^(y)  =  -^'"Z^^^^) 


(95) 
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which  will  have  to  be  discussed  later,  and  signifies  that  the  spectrum 

k~  '   is  actually  in  equilibrium;  however,  for  larger   |  t]  |  the  left  side 
assumes  appreciable  values.   Therefore,  it  is  doubtful  whether  one  will 
obtain  a  sufficient  approximation  if,  in  the  transition  from  (82)  to  (85), 
one  takes  only  the  mean  values  of  the  type  (63)  into  consideration.   How- 
ever, I  did  not  succeed  in  improving  here  the  approximation  or  in 
obtaining  more  than  a  very  crude  estimation  of  (86). 

One  may  perhaps  assiime  for  such  an  estimation  that,  for  large  values 
of  Ti,  the  first  term  of  the  summation  on  the  right  side  of  (9^)  is  much 
larger  than  the  second.   For  in  the  first  one,  the  integral  is  taken  over 
the  region  t]  '  ~  t]  which  probably  contributes  a  great  deal  whereas  in  the 
second,  for  large  t\,    the  factors  /Sg     have  already  strongly  decayed  in 
the  entire  integration  range .   One  may  therefore  attempt  the  assimiption, 
at  least  for  large  t\,    that  the  second  term  of  the  summation  on  the  right 
side  of  (9^)  may  be  neglected.   In  this  approximation  the  time  average  of 
the  first  equation  (86)  then  becomes,  with  use  of  (9^): 


ciTi'g(Ti  +   Tl•)g/Tl•y2/5^  (for     ti  »  0)      (96) 


One  may  utilize  this  equation,  for  instance,  in  such  a  manner  that 
one  assumes  a  plausible  form  for  g('n),  leaving  the  scale  in  the 
Ti-direction  undetermined  at  first,  and  determining  it  subsequently  so 
that  the  equation  (96)  is  valid  as  exactly  as  possible  for  large  t\.      In 
this  way  one  will  describe  the  steepness  of  the  decrease  of  g(Ti)   for 
large   tj  with  some  correctness,  and  precisely  this  steepness  is  decisive 
for  the  value  of  k . 

In  the  practical  execution  of  the  calculation  it  is  expedient  to 
introduce,  in  the  place  of  y  and  f(y),  new  variables 


s  =  y2/3;   cp(s)ds  =  f(y)dy  (97) 

Then  there  applies  as  one  recognizes  from  (87)  (compare  also  (95)): 


cp(i)=  -scp(s)  (98) 
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This  equation  is  based  on  the  fact  that  the  energy  dissipation  from 

the  wave  number  — k  to  the  wave  number  k  coincides  with  the  one  from  k 
a 

to  ak^  except  for  a  factor  qualified  by  the  similitude  transformation. 

Furthermore,  there  then  applies  in  good  approximation 


(p(s)  =  ^(1  -  s5)fl  -  |s5^        (for  0  <  s  <  l)     (99) 


and  for  larger  s  one  can  reduce  cp(s)  by  means  of  (98)  to  the  range 
0  <  s  <  1. 

Figure  4  gives  a  plausible  curve   for  g(ii),  and  in  addition  the  right 

side  of  (96)  as   j   dr]  .  .  .   and  g'(Ti).   The  scale  is  selected  in  such 

h 

a  manner  that  the  two  last  curves  coincide  for  large   i].   Considerable 
differences  then  exist  for  small  r\     but  there  the  equation  (96)  also  can 
no  longer  be  correct.   If  one  substitutes  the  function  g(T])   thus 
obtained  into  (91)  aJ^d.  neglects  the  fluctuations,  there  results  for  k.  : 


K   =  J^  ds(-lg  s)cp(s)     dTig(Ti)g(Tis)  =  0.98       (100) 


This  crude  estimation  therefore  gives  the  correct  order  of  magnitude 
for  K ,   but  the  exact  value  may  well  differ  from  0 . 98  by  as  much  as  a  f ac - 
tor  2.   The  calculations  of  this  section  thus  have  not  led  to  an  exact 
calculation  of  the  constant  k     but  they  did  provide  a  qualitative  mathe- 
matical representation  of  the  processes  on  which  the  energy  dissipation 
is  based.   Perhaps  it  will  be  possible  to  arrive  at  a  rather  exact  experi- 
mental determination  of  k  by  means  of  a  comprehensive  discussion  of  the 
various  experiments  of  Simmons,  Dryden  (cited  before),  Prandtl-'-5,  and 
others  regarding  the  spectrum  and  the  damping  of  turbulence. 

Translated  by  Mary  L.  Mahler 
National  Advisory  Committee 
for  Aeronautics 

5prandtl,  L. :   Proc .  VI.  Intern.  Congr.  f.  Appl.  Mech.  Cambridge, 
Mass.,  1938,  p.  340. 
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Figure  1.-    Representation  of  the  function    w(x). 
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Figure  2.-    The  turbulent  energy  distribution  as  a  function  of  the  wave 

number. 
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Figure  3.-    The  correlation  functions. 
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Figure  4.-    The  correlation  function    g(Ti). 
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